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Structural and Spectral properties of Corona 

Graphs 
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Abstract 

Product graphs have been gainfully used in literature to generate 
mathematical models of complex networks which inherit properties of real 
networks. Realizing the duplication phenomena imbibed in the definition 
of corona product of two graphs, we define Corona graphs. Given a small 
simple connected graph which we call seed graph, Corona graphs are de¬ 
fined by taking corona product of a seed graph iteratively. We show that 
the cumulative degree distribution of Corona graphs decay exponentially 
when the seed graph is regular and cumulative betweenness distribution 
follows power law when seed graph is a clique. We determine explicit for¬ 
mulae of eigenvalues, Laplacian eigenvalues and signless Laplacian eigen¬ 
values of Corona graphs when the seed graph is regular. Computable 
expressions of eigenvalues and signless Laplacian eigenvalues of Corona 
graphs are also obtained when the seed graph is a star graph. 

Keywords: Corona product of graphs, Laplacian spectra, signless Laplacian 
spectra , complex networks , betweenness distribution , cumulative degree dis¬ 
tribution 


1 Introduction 

Network modelling using product graphs is an interesting method to generate 
complex networks which possibly can capture the properties of real world net¬ 
works, for example, Kronecker graphs (see DP)- In (lj, Leskovec et.al. have 
developed both deterministic and stochastic models by using the Kronecker 
product of graphs to generate complex networks which inherit properties of real 
world networks. In [2], the authors had presented generalized graph products 
methodology to generate different types of complex networks. In this paper, we 
propose a model for complex networks by using the concept of corona product 
of graphs. 
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(a) P 3 



Figure 1: Examples of the Corona graphs: (a) A seed graph P 3 (b) Corona 
product of P 3 = G with itself and hence generating G^ (c) The graph of 
another successive corona product of P 3 with previously generated graph GW 
resulting in G^. 


Corona product of two graphs, say G and H , was introduced by Frucht and 
Harary in 1970 [3] is a graph constructed by taking n instances of H and each 
such H gets connected to each node of G, where n is the number of nodes of G. 
Starting with a connected simple graph G, we define Corona graphs which are 
obtained by taking Corona product of G with itself iteratively. In this case, G 
is called the seed graph for the Corona graphs. For instance, the Corona graphs 
generated with G = P 3 is shown in Fig. [T] 

We mention that Corona graphs can be used as models for investigating the 
duplication mechanism of an individual gene as explained for the formation of 
new proteins in [3] and the references therein. A similar phenomenon is being 
observed in Corona graphs but instead of a single node, a unit of seed graph is 
being duplicated and connected to every node of the existing network. Hence, 
Corona graphs can reveal more insights on duplication phenomena and with the 
help of proposed graph spectra, the properties of gene duplication and other 
real world complex networks can be investigated more deeply [2]. 

In this paper, we investigate various structural properties of Corona graphs 
including average degree, sparsity, degree distribution, diameter and cumulative 
degree distribution. We show that 

1. the cumulative degree distribution of a Corona graph decays exponentially 
when the chosen seed graph is regular. 

2. diameter of the Corona graphs is increased by 2 in each iteration of the 
Corona product for any seed graph. 
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Figure 2: Pattern of the adjacency matrices corresponding to (a) G W 
for Fig. 


la 


Dots represent a Vi v . = 1 and white spaces represents a v 


(b) 

= 0 . 


3. cumulative betweenness distribution of a Corona graph follows power law 
when the seed graph is a clique. 

Further, we study the spectra, Laplacian spectra and signless Laplacian 
spectra of Corona graphs obtained by specific seed graphs. We determine 

1. computable formulae for eigenvalues and signless Laplacian eigenvalues of 
Corona graphs generated by a seed graph which is regular or a star graph. 

2. computable expressions of Laplacian eigenvalues of Corona graphs gener¬ 
ated by any seed graph. 

We organize the paper as follows. In Sec. [2] we define the Corona graphs 
and investigate structural properties of Corona graphs. In Sec. [3j we derive 
the spectrum of Corona graphs generated by a regular graph and a star graph. 
Further we determine formula of Laplacian and signless Laplacian eigenvalues of 
Corona graphs generated by a simple connected graph and finally, we conclude 
in Sec. Ill 


2 Corona Graphs 

Let G = ( V ., E) be a graph having the set of nodes V = {iq,..., v n } and E the set 
of edges. The adjacency matrix A(G) = [a ViVj ] of G having dimension |Vj x |F| 
is defined by a VijVj = 1 if ( Vi,Vj ) £ E otherwise a VijV . = 0 when columns and 
rows are labelled by nodes of G. Laplacian matrix and signless Laplacian matrix 
associated with G are given by 

L(G) = D(G ) - A(G) (1) 

and 

Q(G) = D(G) + A(G) (2) 

respectively, where D[G) = diag{d\, d %,..., d n }, di = 'Yhi+j CLij- 

The corona product of the two graphs G\ = (V{G\),E{Gi)) and G 2 = 
(y(G 2 ),E(G 2 )), denoted by G\ o G 2 is obtained by taking an instance of G\ 
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(a) G 

Figure 3: Examples of the Corona graphs: (a) 2-regular seed graph G (b) Corona 
graph for G^. Pattern of the adjacency matrices in (c) G^ (d) . 




and \ V Gl | instances of G 2 and hence connecting the i th node of G\ to every node 
in the i th instance of G 2 for each i [5j. We extend this definition to define Corona 
graphs. Let G^ = G. Given a seed graph G, the Corona graphs generated by 
G are defined by 

G (m+ i) = G (m) o G (3) 

where m(> 0) is a large natural number. For instance, the Corona graphs 
generated by P 3 are shown in Fig |lb| and Fig|Tc| along with the pattern of their 
adjacency matrices in Fig|2a| and Fig [2b] respectively. Similarly, for a 2-regular 
gra ph with 4 nodes, the Corona graph G^ corresponding to G is shown in 
Figj3b| and the pattern of their adjacency matrices of GW and G^ in Fig [3c] 
and Fig|3d| respectively. The following are some observations associated with 
Corona graphs generated by a seed graph of order n. 

1. The number of nodes in G^ m ' > is 

|F (m) | =n(n + l) m . (4) 

2. If \E\ is the number of edges in the seed graph G, then the number of 
edges in G ^ is 

\E^\ = {\E\ + (\E\+n)((n + l) m -l)). (5) 

3. The number of nodes added in i th (i < to) step of the formation of G ^ 
is n 2 (n + l) l_1 . 

4. Connectivity of G t ' rn ' ) : Since G is connected, Corona graphs generated by 
G are connected graphs. Evidently, if seed graph G is disconnected, it 
generates disconnected Corona graphs as shown in Fig[4j This feature is 
not observed in the case of Kronecker graphs given by Leskovec et al. in 
[T| where the authors had taken a multi-graph as a seed graph to ensure 
the connectivity of G^ m \ 
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(a) G 



(b) G« 


Figure 4: A disconnected graph G of 4 nodes having two P 2 graphs as its 
components, along with its 1 st corona product G*- 1 -* resulting in a disconnected 
Corona graph. 


5. Degree sequence of corona graphs: Assume that degree sequence of G ^ = 
G is given by where represents the degree of node 

ii at j th corona product, and x is the total distinct degrees. The degree se¬ 
quence of G is given by {D^\ ... ,n times),..., .. ,n(n + 

l ) m_1 times)}, where D = d , j + inn , = di. + (to — 1 )n,..., 

D [ x) = d, + 1. 

tj 3 


Propostion 1 . Let G= G^ m ~ 1 ' 1 o G be the Corona graph generated by the 
seed graph G. Let \E\ and \V\ = n be the number of edges and nodes in the 
seed graph G. Let \E^ m ^\ and be the number of edges and nodes in G< m ). 

Then, average degree(fk)) of a node in G( m ' > is 


(k) 


2 



Proof. By the definition of average degree of a node, we have 


(k) 


If to » 1, then —> 0. Hence, 

(k) 


2 |£( m )| 

|pM| 

2 ((|£| + (|£;|+n)((n+ !)"-!))) 

n(n + l) m 
2{\E\+n) _ 2 

n (n + l) m 


2 1 + 



□ 
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It follows from above proposition that if the seed graph is a tree such that 
\E\ = (n — 1), then (k) = 2(l + = 2(2 — 1). If the seed graph is a clique 

(K n , having \E\ = Q) edges and n > 3), then (k) ~ (n + 1). 

Density of a network is an important concept to determine the sparsity of 
the network. The density of an undirected network G = (V, E) is defined by [ 6 ] 

d= y (6) 

where \E\ and |U| = n represent the number of edges and nodes respectively in 
G. It is evident that 0 < d < 1. If d <C 1, then the network is called sparse. 
It is noted in [7j that many real world networks have d £ [10~ 5 ,10 -1 ]. In [ 8 j, 
the authors remarked by citing the example of neural networks that sparsity 
saves the energy of a network without effecting its functionality. The network 
sparsity is also observed in other biological networks like metabolic network (0 
and the references therein). 

Propostion 2. Let \E\ and n be the number of edges and nodes respectively 
in seed graph G. Let \E^ m ^\ and be the number of edges and nodes 

respectively in G^") = G l ' m ~ 1 ' > o G. Then, for m 1, the Corona graphs 
G< "b are sparse as d <C 1 . 

Proof. Using equation ([b]), the density of G is given by 

- (Vfi) 

2 (|£’| +n) n 

~ n(n(n + l) m - 1 ) ("("+i) m )' 

If to 1, then —> 0 and also 1. Hence, the result 

follows. □ 

2.1 Degree Distribution 

Now, we consider the degree distribution of Corona graphs. In a network G, 
let P(k) denote the fraction of vertices having degree k and hence the degree 
distribution is the probability distribution of the different degrees in the whole 
network mm ■ The total instances for a particular degree k in G^ is given 
by 

N k = n (i) + n ^ fc £) ( 2 ) + n {n+l)S k D (.3) + ■ ■ ■ + n(n+ l) (m (^) 

■ , V ’ 'i ’ '3 ' 'i ’ '} J 
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Figure 5: Degree Distribution for seed graphs (a) K$ with and G^ having 

49152 and 12288 nodes respectively, (b) P 3 with G^ and G^ having 49152 
and 12288 nodes respectively. 


where 5 k ^ is the Kronecker delta function. Therefore, the degree distribution 

’ % 3 

P(k) of Corona graphs is given by 


m = 


E n 

.7 = 1 


S k, D™ 


nS kM 2) 


n(n ■ 


1 )' 5 fc.D (3) 


+ ... + n(n + 1)^ 




nln - 


i) r 


( 8 ) 


The degree distribution for K 3 and P 3 are shown in Figj5] A crucial obser¬ 
vation from the degree distribution is that there is huge multiplicities of degrees 
in the Corona graphs and this can also be confirmed from the Fig|5a[ Fig|5b| 
It is similar to the observation as seen for Kronecker graphs Tj. The figure also 
shows that the degree distribution follows similar type of curves for different 
instances of Corona graph i.e. for for each i £ [1, m\. The figure also shows 
the fat tailed degree distribution in both sub-figures of Fig(5] The fat tailed 
distributions are found abundantly in real world networks like data traffic on 
internet, return on financial markets etc. as mentioned in | 12 j and the references 
therein mm ■ It should be noted that the degree sequence of a Corona graph 
comprises of discrete values, so the degree distribution of a Corona graph could 
be properly investigated with cumulative degree distribution. 


P c (k) = P(k') 

k’>k 


which is the probability that the degree is greater than or equal to k. Here, 
P(k) and P c (k) are the discrete and continuous degree distribution respectively. 
Thus, we have the following result. 

Theorem 1. Consider the Corona graph G^ = o G generated by the 

seed graph G which is r-regular. Then, the cumulative degree distribution of 
Corona graphs decay exponentially. 
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Figure 6: Cumulative degree distribution of G^ generated with the seed 
graph = K 3 where k is degree of nodes and P c (k ) is cumulative 
probability at each k. 

Proof. The generating methodology of Corona graphs for t > 1 implies 

ki(t) = ki(t- 1) +n (9) 

where kjfi) is the degree of the node i in G^\ Since the degree of any node is 
r + 1 while it gets added during formation of , it implies that 


hit) = r + 1 + n(t - U). 


( 10 ) 


Now, cumulative degree distribution can be obtained using equation ( 10) as 
follows 

P c (k) = Y P(k ) = P(t < t = t H-—- 

< ^ 71 


k >k 


Hence, 


P c (k) = 


Since t = t + r+ * k , hence 
Pc(k) 

Since k > r + 1, the result follows. 


V 

|F(*)| 

t'— 0 

|V(‘')| 

\y(t)\ 1 |y(t)| 


E 


1 , ((n + l) T -l) 


[n + l) 4 (n + l) 4 

(n + l) r_t 


= (n + 1) 1 


( 11 ) 

□ 














The Figgis showing cumulative degree distribution of G^ with seed graph 
G(°) = A 3 . The exponential decay of degree distribution is also observed in m 
when the derived resultant graph has same number of offsprings for all nodes. 
The proposed model in m has also exponential degree distribution. 

2.2 Diameter 

Diameter of a graph is the longest shortest path in the graph. We show in the 
following theorem that as a Corona graph grows, the diameter also increases. 

Theorem 2. Consider the Corona graph G^ m ^ generated by a seed graph G 1 ' 0 ' 1 
with the diameter . Then, the diameter of Gis D + 2m. 

Proof. We prove it by induction method. Let i, j be the end nodes of the 
diameter D in G^°\ In G 1 = G ^ o G^°\ a single instance of G will be 
attached to every node(according to definition of corona product) including 
i and j and also each node will be connected to each node of G 1 ' 0 ' 1 . Now, 
the diameter(.DW) of the graph will be elongated as D + (an edge from i) + 
(an edge from j) = ZA°) + 2 . 

For inductive hypothesis, let the diameter of the Corona graph G^” 1-1 ) be 
_ £)(o) 2 ( to _ x) and ( i ,j ) are the two end nodes of Again, 

in G^ m \ a single instance of G^ will be attached to every node of G^ m_1 ^ 
including (i ,j ) and hence D^ m ' 1 in G (m) 

is +2(m— l) + (an edge from i ) + 
(an edge from j ) = D + 2 to. □ 

2.3 Betweenness Distribution 

Betweenness evaluates the importance of a person in a social network (see ITTl 
and the references therein). It is also being assessed in the models studying the 
cascading failures as the load on the i th node |T0] , p~8] , [19 j . Mathematically, it 
is formulated as 

* = £ —21 (12) 

tidl, ”><• 

where hi is the betweenness centrality of the i th node, <Jjk(i) is the number of 
the shortest paths between nodes j and k through i, and Ujk is the total shortest 
paths between j and k. 

It is difficult to find out all the pairs of shortest paths for crjk(i ) and ajk 
for the Corona graphs generated by a seed graph. We observe that there is a 
distinct shortest path in a Corona graph between all pair of nodes when the 
seed graph is a clique. There are two reasons of this observation-(l) each node 
in a clique is connected to each other node of the clique and hence making a 
unique shortest path between each other. ( 2 ) each seed graph G ^ is connected 
in the i th corona product to any of the node rij of G*- 1-1 ) and this node nj is 
connected with a unique path to all nodes of this G^ as well as it is acting as 
a unique bridge (and hence a unique path) between the nodes of G^ and the 
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rest of the nodes of G (* ^ and G^\ Hence, when G^ = K n (where n > 2), 
eqn.(13) can be modified as 



(13) 


j^i^k 


Now, for assessing the number of shortest paths, we use reasoning similar to 
mm- Let all the nodes attached to node i be termed as offsprings. Similarly, 
the nodes attached to the offsprings as well as the nodes attached recursively 
to all the offsprings in the subsequent corona products G^ (corresponding to 
node i) are being termed as offsprings. 

Theorem 3. The betweenness distribution of a Corona graph generated by a 
clique as the seed graph follows power law with exponent approximately 2. 

Proof. There are two types of shortest paths in the Corona graph G^ generated 
by the seed graph K n through node i: 

1. A = Total shortest paths between offsprings to the rest of the nodes of 
the graph. 

2. B = Total shortest paths between offsprings. 

such that b\(r) = A + B where b\(r) represents the total number of shortest 
paths through the node i which was added in r th (r < t) step during the 
formation of G^\ For f > 1, 


6t(r)«>l=|yW|(|yW|-|yW|) 


where |U( T )| and |V^| represent the total offsprings in G (where r > 1 ) 
(corresponding to node i ) and number of nodes in G W such that 


T 


|W T )| = ^n(n + l y - 1 


= ((n + iy-1). 


(14) 


and 


|T (t) | =n(n+l)‘. 


(15) 


Using equations ( 141 and ( 15), we get 


bt(r) = ((n + 1 ) T - l)(n(n + l ) 4 - n + l) r + 1 ). 

For t 1 and r > 1, we get 

b t {r) » n{n + l) t+T ~\ (If 

Let cumulative betweenness distribution be denoted by P c (b). Then, 


(16) 
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Figure 7: Cumulative power-law betweenness distribution of G ^ generated 
with seed graph G ^ = K 3 as an instance of a clique. The betweenness values 
are generated with igraph [22]. 


Thus 


Pc(b) 


t-T -\-1 

E 


M= 1 


\v^\ 

|yw| 


n(n + 1)* r+1 
n(n + l) 4 
(n + l) 1-r 

|yW] 

n(n + l) t+r_1 
(n(n + l) t+T_1 ) _1 . 


(n(n + l) t+T-1 ) 1_7l> = (n(n + i)*+ T - 1 )- 1 


which implies 

7b « 2. (17) 

□ 

The Fig. [T] is showing cumulative power-law betweenness distribution for 
G generated by the clique = K 3 as the seed graph. The betweenness 
distribution featuring power-law was also observed in [23], [24], [25]. In [ 2D] , 
the authors analysed a family of deterministic recursive trees having exponen¬ 
tial decay in cumulative degree distribution as well as cumulative power-law 
betweenness distribution. 
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3 Spectra of Corona Graphs 

Let Gbe the Corona graph generated by the seed graph G. The adjacency 
matrix A(G^) associated with G ("b is given by 

A{G^~V) l^ (n+1) »-i®/n (n+ l)»-i' 

l-n(n+ l) rn ~ 1 & In(n+ 1)™ -1 ^-(^) ® ■Aifn+l)™' -1 


A(G (m) ) = 


where ^4(G^ m 1 - ) ) is the adjacency matrix of G^ m 1 \ /„( n +i)m-i is the identity 
matrix of order n(n+ l ) 1 ™^ 1 and l n ( ra +i)™-i is the all-one vector of order n(n + 
l) m— 1 . We denote spectra of the corona graph G("b i.e. spectrum of A(G < ' m ' > ) 
by 

a(G^) = { Ai,A 2 ,...,A n(n+1)m } (18) 

where Ai < A 2 < ... < A„( n+1 ) m and spectral radius of A(G < ' m is denoted as 
p(G( m )). The spectrum of corona product of two graphs G = Gi o G 2 where G\ 
is any graph and G 2 is a regular graph, and the Laplacian spectrum of corona 
product of any two graphs are provided by Barik et. al. in [5]. Inspired by their 
work, we derive cr(G' m )) when the seed graph G is regular. 

In the next theorem, we provide a(G( m ' > ) in terms of the eigenvalues of the 
seed graph = G. 

Theorem 4. Let G^ = G be a regular graph such that a(G) = (/xi, P 2 ,... ,p n = 
r ) (where, pi < p ,2 < ... < p n = r) and spectral radius of G be p{G). Then, 
er(GM) i s given by 


(a) 


hi + r E^o 1 2 ° ± (e” - ! 1 - + \j ((?’ — hi) ± El^i 1 z c) 2 + 2 2; 

A* = - V - 

2 m 

fj(G( m )), with multiplicity 1 for i = 1 ,..., n(n + l)" 1-1 
where, 


Zk -1 = Zk —2 + \/(( r - ± Zfe- 2) 2 + n.2 2 ( fc_1 ), 2 < k < m, z 0 = 0. 

(b) /ii £ a(Gwith multiplicity n(n + l) m_1 /or i = 1 ,..., n — 1. 
The spectral radius of G^ m ' is given by 



£ 


p(G (m) ) = 


+^e:= 0 1 2“ 


+ (E c =i 


+ ^/((r-^)-E?Ji z c) 2 + 2 2m -n) 


where z c is defined above. 

Proof. Let 1„ be the all-one vector of dimension n. Let Z\,Z^,...,Z n be 
the eigenvectors associated with eigenvalues Ai, A 2 ,..., A n of A(G) respectively. 
Then, the spectrum of G*- 1 - 1 = G^ o G^ (as given in Theorem 3.1 of [5] ) is as 
follows 
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(i) tJ,l+r± V^ Vi) + 4 ” e fr(G'^) with multiplicity 1 for i = 1 Let A 

and Xi be their representations. Then, eigenvectors corresponding to them 
are 

Zi \ f Zi 

Xyry(l n ® Zi)J 7 y^._ r (ln ® ^») 

(ii) fa £ a(G^) with multiplicity n for i = 1 ,..., n— 1 and the corresponding 
eigenvector is 

( ° 

\^Zi G ej 

where j £ [1 , n] and ei,e 2 ,... ,e n are the unit vectors of standard basis 
for n x n. 


Now, the eigenvalues for G= G^ m ^ o G are 

Mi + r EZ~o 2“ ± (e^i 1 ^ + /((r-M ! )±Eri 1 ^) 2 +2^.n) 


(a) A = 


a(G < ' m ' ) ), with multiplicity 1 for ? = 1,... ,n(n + l )” 1 1 corresponding to 
the eigenvectors 


C// m) = 


t/, 


(m-l) 


, A. -r \ 






where, 


= (>-« dr- 2 ') Evb (i,• j 


l a £ [1,2 n{n + l) a - 1 ] such that A i a is an eigenvalue of G^ a \ Zk -i = 
Zk -2 + V ((V - Mi) ± Zfc- 2) 2 + n.2 2 ( fc ” 1 \ 2 < k < to, z 0 = 0. 

(b) /ii £ a(G^), with multiplicity n(n + l) m_1 for i = 1 ,..., n — 1 and the 
eigenvector corresponding to these eigenvalues is 

( ° ) 

\Zi ® ej) 

where j £ [l,n(n + l) m_1 ] and e ? - are the unit vectors of standard basis 
for n(n + l) m_1 x n(n + l) m_1 . 

□ 


Fig. [ 8 ] is showing the eigenvalue distribution of Corona graphs, with the seed 
graph = x 3 _ it i§ the ratio between the multiplicity of eigenvalues and the 
number of nodes. We derive the spectra of star graphs Sk with k > 3 which is 
a special case of irregular graph. 
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O • n=12 

O n=48 1 
n=192 


.0 

- 

0 : 

0 \ 


/ \: 6 



j y 




(a) 


(b) 


Figure 8: Eigenvalues distribution with the seed graph = K 3 for (a) 

G^,G^ and having 12,48 and 192 nodes respectively, (b) having 
3072 nodes. Here, P m (A) is showing the probability distribution for G( m \ 


Theorem 5. Let k > 3 be an integer. The spectrum of the graph Sk°Sk consists 
of the following eigenvalues 


(a) X z = Xf + ^ £ cr(G(■*■)) with multiplicity 1, and 

(b) 0 £ a(G^) with multiplicity k(k — 2). 


where, \ z are 

8 27T+8 4-rr+e 

3’ 3 ’ 3 


the eigenvalues for G ^ with z = 
as shown in following sub-expressions 


1,2,3 for the 3 angles i.e. 


X z = | cos 


yn+8 /,,2 I tap o\ a _ m „-l ( 2/*?+/i<(18-9fc)+(54fc-54) \ 

-H-VK+W 3),0-cos ^ 2(/i?+(6fc _ 3))i j F 

where y = 0,2,4. Here, A z £ 


2 Mi 2y'(6fc-3) 2 1 /(6fc-3) 

3 3 j Pi t 3 


Proof. Let 1„ be the all-one vector of dimension n. Let Z 3 , Z 2 , ■ ■ ■, Z n be the 
eigenvectors corresponding to seed graph’s eigenvalues p\, P 2 ,..., p n respec¬ 
tively. The eigenvalues corresponding to G^ are X z = X[ + ^ £ a(G^ l ' > ) for 
z = 1, 2, 3 and the eigenvectors corresponding to them are as 


Ai + 1 

Xf-k+l 


Z t 

(lr 


8 > Zi) 


GW has 0 as its one of the eigenvalues and its multiplicity is of k(k — 2). □ 

Corollary 1. Let G be the seed star graph S k for each k > 3 such that a(G) = 
{pi, P 2 , ■ • •, p n }- Let m > l.Then cr(G^ m ^) is given by 

(a) X zA =*? + $,..., X ZjTn = Er=o W) jX m-i + (|r _1 (^) G a(GW) 
with multiplicity 1 for each of them, 
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(b) 0 € a(G^ m ^) with multiplicity k{k — 2 ){k + l)l m 1 ^. 

where, X z j are the eigenvalues for G^ m ^ such that j represents j th corona prod¬ 
uct with z = 1, 2, 3 for the 3 angles i.e. ; 47r + ( b as shown in following 

sub-expressions 

xf = § co S y-^^{Y!- = \{\) J - lx t-j + (D'-VO 2 + (6 k 3), 


| cos 

9 m = cos -1 
yn 


yn+di 


Vvf + ( 6 k - 3) 


U^)) 2 +(6fc-3))2 


2 ((Er= _ o 

where m is the m th corona product and y = 0 , 2 ,4 for the three angles. 

2m» __ 2j v /(6fc-3) 2jy/(6k-3)~ 


Here, X z j £ 


-Vi 


3J 


i Vi 


, where j = 1 ,..., m. 


Proof. Let 1„ be the all 1 vector of dimension n. The proof follows by using 
similar arguments given in the proof of Theorem [4] and Theorem [5] The 
eigenvectors corresponding to the eigenvalues X z j, z = 1,2,3, j £ [1, m] are 
given by 

[/('-!) \ 


u (l) = 


/ -\~k —1 



. K ,• -fc+i 


C/('- 2 )) 







A 2 i j, are the eigenvalues for G^ l \ 1 < l < m. 

G^ has 0 as its one of the eigenvalue and its multiplicity is of k(k — 2 )(fc + 

□ 


In Fig. [9j we plot the distribution of eigenvalues generated with the seed 
graph G^ = S 3 . 

The Laplacian matrix L{G^ m ' > ) associated with G( m ' > has the form 


L(G (m) ) 


L(G l ' m 1 ' ) ) + nl„( rt+ 1 )m-l —® Infn+l ) 1 "- 1 
. Infn+l) 7 ™ -1 ^ In(n-f-l) m_1 (^(G) T ® In(n+l) m_1 . 


where L(G^ m is the Laplacian matrix of G^ m 1 \ I n and I n / n are the 

identity matrices. We denote the Laplacian spectra L(G^ m ' > ) of G^ by 

S(G^) = {A,, A 2 ,..., A n(n+1) (m)} (19) 


where 0 = Ai < A 2 < ... < A n ( n+1 vm). In the following theorem, we determine 
the elements of S(G^ m ' ) ) in terms of the Laplacian eigenvalues of the seed graph, 
G ( °) = G where S(G) = {0 = vx,u 2 , ..., u n } (where, v\ < v 2 < ■. ■ < v n ) is the 
Laplacian spectra of G. The algebraic connectivity of a graph is defined as the 
second smallest eigenvalue of L(G) [26]. 
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Figure 9: Eigenvalues distribution with the seed graph G^ = S 3 for (a) 
GC 1 ) ,G f ( 2 ) and G^ having 12,48 and 192 nodes respectively, (b) G^ having 
3072 nodes. Here, P m (A) is showing the probability distribution for G( m \ 


Theorem 6. Let G be a simple connected graph. We denote the algebraic 
connectivity of G and G( m ' > by a(y 2 ) and a( A 2 ) respectively. Then, Laplacian 
spectra S ( G ( m ' > ) of G is given by 


(a) 


^+(n+i)z?=- 0 1 r±i:Ziz i 

2 m 

1 ,..., n(n + l ) m_1 . where, 
z 1 = \/(i'i+n + l ) 2 - 4i/, , 


G 


5 (G( m )) with multiplicity 1 for i = 


Zk 


i 


(y% + (n +1) 


k -1 


k -1 


± ^ Z i) 2 


2=0 


for 2 < k < m. 


k—2 k -1 

2(fc+ 1 )(j / . + ( n + 1) E 2i± E*<> 

i=0 i—1 


(b) i/j+le S(G < ' rn ' ) ) with multiplicity n{n + l) m 1 /or i = 2,... ,n. 
Hence, the algebraic connectivity of S(G^ m ' ) ) is 


a(A 2 ) 


^ + (n + l)E”o 1 2 t 
2 m 


E ra 
i —1 ■ 


< 1 


where Zi can be defined as above. 


Proof. Let 1„ be the all-one vector of dimension n. Let X\, X 2 ,..., X n be the 
eigenvectors of corresponding to the eigenvalues v\, z/ 2 ,..., n«„ respectively of 
L(G). Then, the spectrum of G W = o G® (as defined in Theorem 3.2 of 
0 ) is as follows 

(i) I ' i+n+ 1 ± VG^+ n + 1 ) 4t/ » g S^G^) with multiplicity 1 for i = 1,..., n. Let 

/ii and fii be their representation. Then, eigenvector corresponding to 
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them is as 


Xi \ ( Xi \ 

X~“(lra <8 Xi) J ’ y \-( Li (!«■ ® X{) J 

(ii) Vi +1 £ S(G ^) with multiplicity n for * = 2 ,..., n and the corresponding 
eigenvector is 

( 0 

\X i <g> ej 

where j G [l,n] and e\, ..., e n are the unit vectors of standard basis 

for n x n. 


Now, the laplacian eigenvalues for G^ m ' ) 

^+(n+i)E^ 1 2‘±E;: 1 ^ p 

v / 2 m ^ 

corresponding to the eigenvectors 


= oG is 

S(G^ m ^) with multiplicity 1 for i = 


where, 




U, 


(m-l) 


lm-1 


<U> 


i-i )' 


Min-l) 



u l™~ 2 2) 


c/ (1) - 

J U U — 



l a G [1, 2?i(?r + 1)° x ] such that Aj a is an eigenvalue of G^ a \ 
zi = ij{vi+n + l) 2 - 4 1 /, , 


— 


\ 


m— 1 m—1 


m—2 m— 1 


K + (n + 1) ^ 2* ± ^ ^) 2 - 2(™+!) + (n + 1) ^ 2* ± ^ Zi ) 


0 2=1 


2—0 2—1 


for m > 2 . 

(b) i/j + 1 £ 5(G^ m ^), with multiplicity n[n + l ) m_1 for * = 2 ,, n and the 
eigenvector corresponding to these eigenvalues is 


X, 


where j G [l,n(n + l) m x ] and e :l are the unit vectors of standard basis 
for n(n + l ) m_1 x n(n + 1 ) TO_1 . 


□ 
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(a) (b) 

Figure 10: Laplacian eigenvalues distribution with the seed graph G^ = I \ 3 for 
(a) G^jG^ and G® having 12,48 and 192 nodes respectively, (b) G^ having 
3072 nodes. Here, P m (A) is showing the probability distribution for G( m \ 


Fig. [TO] is showing the laplacian eigenvalue distribution, with the seed graph 
G ( °) = K 3 . 

The signless Laplacian matrix SQ(G^ m ^) of Gis of the form 

g ^Q(m)) _ Sq(G^ ^) + nl n ( n+ i)m-l l„( n+1 )m-l <S) Infn+l)™- 1 

. ln(n+l) m_1 ® In(n+l) m_1 (^q(^) 4” ^n) & In(n4-l) m_1 . 

where Sg(G ( - m ~ 1 ^) is the signless Laplacian matrix of G^ m_1 \ I n and 7„( n +i)m-i 
are the identity matrices. The spectrum of signless Laplacian of <5(G^ m - ) ) for 
G (m ) is denoted by 


Sq(G^) = {A 1; A 2 ,..., A n(n+1) ( m) } (20) 


where Ai < A 2 < ... < A ra ( n+ i)m. Recently, there is a lot of work on signless 
Laplacian matrices and their Q —spectra as in E2,EBj,E2] and their authors 
think that it is more useful and hence a lot of work is going on to find their 
usefulness as mm- We will define the Sq (G^) inspired by the work of 
m on two graphs and here taking the seed graph G as the r-regular graph. 
In the following theorem, we derive the elements of Sq(G^) in terms of the 
signless Laplacian eigenvalues of a regular seed graph, G^ = G such that 
Sq{G) = (q u q 2 ,..., q„ = 2r) (where, <?i < q-i < ■ • ■ < q n = 2 r). 

Theorem 7. Let G be a simple connected graph. Then, Sq (GM) is given by 


(a) 


x _^ + n EZo 1 * + r Er=i 2 l + EEo 1 2‘ ± E£=i 

Aj — -—- 

L o m 

of 1, for i = 1 ,..., n(n + l) m_1 
where, 


with multiplicity 
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3 - 1 3 — 1 3 ~ 2 i - 1 

(ft + n^2 i + r(^ 2* - 23) + (^2 2 i - 23- 1 ) ± ^ft ) 2 + 2 2 Ln 
2=0 2=1 2—0 2—1 

/or j = 2,..., to and z\ = \J ((ft + n) — (2r + l )) 2 + 4n. 

(b) qj + 1 with the multiplicity of n{n + l ) m_1 /or / = 1 ,..., n — 1 

Hence, spectral radius of Sq ( G ( m ' > ) is 



q(S Q (G™)) = 


ft+«Er=o 1 2 i +^r=i2 i 


E m—1 C)i 
2 — 0 ^ 


E m 

3=1 Z 3 


where Zj is defined as above. 


Proof. We can prove the Part (a) of the theorem by induction and hence, the 
q(SQ(G ( ' m ' ) )) will be followed from the proof which is as follows 
Base case: For j = 1, G^ = G^ o G , the Sq(G can be defined according 
to Theorem 3.2 of m as 

(i) ^+"+ 2 ^i±v/((^+»)-(^+D) 2 +^ € Sq{g( i) } with multiplicity i for i = 


(ii) ft + 1 € Sq(G'- 1 >) with multiplicity n for i = 1 ,..., n — 1 

Inductive hypothesis: Let j = m — 1, G^” 1-1 ) = G^" 1-2 ) o G, the S^G^" 1-1 )) 
can be defined as 


v—^vm—2 

(a) A. = + 

ity of 1 , for i = 1 ,. 
where, 


■rES 1 ^ 


E m—2 
2 — 0 


2 m —1 

.., n(ra + l ) m_2 


2 * ± W m_1 ft 

- 2 ~^~ 1 with multiplic- 


1- i j-i j- 2 j-i 

{q i +nJ22 i + r(^ 2 l - 2 ■?) + 2 * - 2 J- 1 ) ± ^ ft ) 2 + 2 2 Gn 

2 — 0 2=1 2=0 2=1 

for j = 2,... ,m — 1 and z\ = \/((ft + n) — (2r + l )) 2 + An. 

(b) ft + 1 with the multiplicity of n{n + l ) m-2 for i = 1 ,..., n — 1 

Inductive step: For j = m, G^ m ^ = G^ m_1 ^ o G, the Sq(G^) can be obtained 
by substituting the eigenvalues of SQ(G^ m ~ 1 ^) in place of ft of the base step 
and we will get the eigenvalues as stated in theorem. 

The other eigenvalues are ft + 1 with the multiplicity of n(n + l ) m_1 for 
j = 1 ,... ,n — 1 . □ 

Fig EH is showing the signless laplacian eigenvalue distribution, with the 
seed graph G= K 3 . 

Consider star graph Sk which is an irregular graph. In the theorem below, 
we determine explicit formula of signless Laplacian elements of G ( m ' > = . 
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Figure 11: Signless laplacian eigenvalues distribution with the seed graph G^ = 
K 3 for (a) G^\G^ and G^ having 12,48 and 192 nodes respectively, (b) G^ 
having 3072 nodes. Here, P m ( A) is showing the probability distribution for 
G( m ). 


Theorem 8. Let Sq(G) = {qi,q 2 , ■ ■ ■ ,q n } with <Zi < 92 < ■ • ■ < q n - Then, 
Sq(GW) is given by 

(a) X z = Xf + 9 ^+ 2 k +2 g Sq(G^) with multiplicity 1, and 

(b) qj + 1 £ Sq(G^) where for j = 2 ,... , n — 1 with multiplicity of each of 
them as k 


where, \ z are the eigenvalues for G ^ with z = 1,2,3 for the 3 angles i.e. 

2?r 3 l ~ e , 47 y~ e as shown in following sub-expressions 
X z = | cos yjq? + qi(k - 2) + (k + l) 2 , 


Q _ -i( 2 9 ^ + (3fe —6)q, 2 —3(fe 2 —fe—2) 9i + (T0fc —94—12 Ea7i ( a + 2 )(fc-a-l) \ 

(7 — COS I - 3 - I — 

2{qf+qi(k-2)+{k+l) 2 )2 


V 7T 


where y = 0,2,4. Here, A z e 


, M . y-l^+v 7 ^-^ 2 ) , Y+A+s/IB-A 2 

m I o 1 O T o 


w/iere H = (/c — 2), H = (fc + l) 2 , Y = (2k + 2). 


Proof. Let l n be the all 1 vector of dimension n x 1. Let Y\,...,Y n be the 
eigenvectors corresponding to seed graph’s eigenvalues q\,,q n . The eigen¬ 
values for G W are A- = Xf + gi+2fc+2 £ Sq(G^) for 2 = 1,2,3 and the 
eigenvectors corresponding to them are as 


Yi 

qi-k+l f 1 

g 2 - 9 i (fe+ 2 ) + (fc+l)l- L ^ 




1 j + 1 e S'q(G^) are the other eigenvalues for j = 2 ,..., n— 1 with multiplicity 
of each of them as /c. □ 
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Corollary 2. Let Sq(G) for star graph Sk for each k > 3 is {q±, < 72 , • ■ •, q n } 
with qi < q 2 < . ■. < q n . Then, Sq(G is given by 


qi+2k+2 \ _ 1/ I s ) 

3 ’ 

£ SQ(G < ' m ^) with multiplicity l, and 


(a) A z ,i — Xf + ——g- 1 • • • j A Zi77 j — y ‘j—n ( 3 ) J Xm-j + ( 3 ) m Qi + (2fc + 


(b) qi + 1 £ Sq(G ( m ' ) ) where for i == 2,... ,n—l with multiplicity of each of 
them as k(k + . 

where, X z j are the eigenvalues for G^ (for j th corona product) with z = 1,2,3 
for the 3 angles i.e. y, 27T + 0J ) 47r + 8j as shown in following sub-expressions 

Xf = ^/(A) 2 + A + (k + l ) 2 _ 

A"f = | cos y7r 3 +e \/— 2 ) + (k + l) 2 , 

where A = + (iT^qi + (2k + 2) Z?Ji(l) j ) 

2A^ m _ 1 +(3fe-6)A^ m _i-3(fc 2 -;s-2)A a , m _i + (70fe-94-12X;^£?(o+2)(fe-a-l) 


2(A 2 m _ 1 +A., m _i(fe-2) + (fc+l) 2 ) 2 


0 = COS 
J/7T 

where y = 0,2,4 and A Zim _i is as defined in part(a) of corollary. Here, 


A z 1 £ 


, 2gi , y-(A+V4B-A2) Y+A+n/4B-A2 

Hi ' 3 ' 3 1 HI ' 3 

^(-2 E”Li py + E™ 1 3"*), qi + m( 

where A = (k — 2), B = (k + l) 2 , Y = (2k + 2). 


3 m 


; • ■ ■ > E 


—qi + 


Y+A+V4B-A 2 


Proof. Let l n be the all 1 vector of dimension nxl. Let Yi,..., Y n be the eigen¬ 
vectors corresponding to seed graph’s signless laplacian eigenvalues qi,... ,q n . 
Before presenting the eigenvectors corresponding to the part (a) of the corol¬ 
lary, let be the eigenvector corresponding to (where i £ [ 1 , m]) and 
li £ [l,3n(n + 1) I_1 ] such that A = g 2 _ g .(E^)+(fc+i) s ig n ^ ess laplacian 

eigenvalue (defined by above expression) corresponding to P/.E The eigenvec¬ 
tors are as follows 

p(m-l) 

Mm -1 

A (l„®p/ m m - 1} ) 

where 



p(m-l) 

Mm-1 


P 


(m— 2) 


lm -2 


, A ( 1 ,; 


3 (m-2)> 




E \ 

(In ® E)y' 


q.i + 1 £ o^Gl 1 )) where for i = 2 ,... ,n—l with multiplicity of each of them 
as /c(/c + l)( m_1 ). □ 

In Fig. [l2j we plot the distribution of signless Laplacian eigenvalues gener¬ 
ated with the seed graph G^ = S 3 . 
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(a) 


(b) 


Figure 12: Signless laplacian eigenvalues distribution with the seed graph G ^ = 
S 3 for (a) gC 1 ) .gC 2 ) and G ^ having 12,48 and 192 nodes respectively, (b) G^ 
having 3072 nodes. Here, P m ( A) is showing the probability distribution for 
CH. 


4 Conclusion 

We proposed a model for generation of complex networks inspired by the phe¬ 
nomena of duplication of genes. We defined Corona graphs by taking corona 
product of a simple graph, which we call a seed graph, finite number of times. 
We determined structural properties of the Corona graphs including cumulative 
degree distribution for any seed graph and cumulative betweenness distribution 
when the seed graph is a clique. We determined spectra, Laplacian spectra and 
signless Laplacian spectra for corona graphs when the seed graph is regular. We 
also derived the spectra and signless Laplacian spectra of corona graphs when 
the seed graph is a star graph. 
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